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ABSTRACT 

A next generation reusable launch vehicle (RLV) 
will require thermally efficient and light-weight 
cryogenic propellant tank structures. Since these tanks 
will be weight-critical, analytical tools must be 
developed to aid in sizing the thickness of insulation 
layers and structural geometry for optimal performance. 
Finite element method (FEM) models of the tank and 
insulation layers were created to analyze the thermal 
performance of the cryogenic insulation layer and 
thermal protection system (TPS) of the tanks. The 
thermal conditions of ground-hold and re-entry/soak- 
through for a typical RLV mission were used in the 
thermal sizing study. A general-purpose nonlinear 
FEM analysis code, capable of using temperature and 
pressure dependent material properties, was used as the 
thermal analysis code. Mechanical loads from ground 
handling and proof-pressure testing were used to size 
the structural geometry of an aluminum cryogenic tank 
wall. Nonlinear deterministic optimization and 
reliability optimization techniques were the analytical 
tools used to size the geometry of the isogrid stiffeners 
and thickness of the skin. The results from the sizing 
study indicate that a commercial FEM code can be used 
for thermal analyses to size the insulation thicknesses 
where the temperature and pressure were varied. The 
results from the structural sizing study show that using 
combined deterministic and reliability optimization 
techniques can obtain alternate and lighter designs than 
the designs obtained from deterministic optimization 
methods alone. 

INTRODUCTION 

The development and fabrication of reusable 
cryogenic tanks is one of the significant technical 
challenges in the development of an operable RLV [ 1, 


2]. Large (30-ft. dia.) reusable cryogenic tanks will be 
required to contain the LH 2 and LOx for an RLV. The 
development of flight-weight reusable cryogenic 
propellant tanks is critical for an RLV because the 
tanks, as shown in Figure 1 for a generic RLV, may 
comprise as much as 35 percent of the dry structural 
weight of the vehicle and occupy 70 percent of the 
internal volume of the vehicle. The cryogenic tanks of 
an RLV must not only function as pressure vessels at 
cryogenic temperatures, but the tanks also must carry 
primary structural loads and support the thermal 
protection system (TPS). An example of a section of 
an integrated tank system design is shown in Figure 2. 
In this concept, an aluminum 2219-T87 (A1 2219) wall 
with cryogenic insulation (Cl) serves as the pressure 
vessel, primary structure, and TPS support structure. 
Analytical tools must be developed to size the through- 
the-thickness insulation and structural geometry of the 
tank structures. With these tools, various combinations 
of insulation concepts and tank structures could be 
efficiently analyzed and optimized in thermal and 
structural trade studies to develop robust and light¬ 
weight integrated structures for an RLV. 



surface 

Figure 1. A generic RLV with an aft-located LH 2 
tank. 
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The objective of this paper is to present the 
analytical tools that were developed and results for the 
thermal sizing of insulation for the Cl and TPS, and 
structural sizing of the geometry of the isogrid 
stiffeners and skin thickness for an aluminum cryogenic 
propellant tank of an RLV. A nonlinear finite element 
method (FEM) analysis program capable of performing 
thermal analyses using temperature and pressure 
dependent material properties was utilized to analyze 
combinations of candidate Cl and TPS. Nonlinear 
structural analysis, deterministic optimization 
techniques, and reliability optimization techniques were 
used to size the isogrid stiffening geometry and skin 
thickness for an aluminum cryogenic tank. The 
analytical tools, methods, and results of the thermal 
sizing study will be presented first, followed by the 
analytical tools, methods, and results of the structural 
sizing study. At the end of the paper, the results from 
both sets of analyses will be combined. 



Figure 2. An example of an RLV metallic 
integrated tank concept with an aluminum isogrid 
tank wall, foam-filled honeycomb insulation, and 
conformable reusable insulation (CRI) TPS. 

THERMAL SIZING STUDY 

A thermal sizing study of two integrated tank 
design concepts was conducted to develop the 
insulation thickness of the Cl and TPS and to 
maximize the thermal efficiency for two load cases: 1) 
ground-hold and 2) re-entry/soak-through. Nonlinear 
transient heat transfer analyses for the two cases were 
performed using ABAQUS, a nonlinear FEM analysis 
code. 

Approach 

The criteria used for sizing the thickness of the Cl 
and TPS for an LEE tank for the ground-hold case were 
to: maintain propellant quality (average heat flux 
through-the-thickness less than 45 BTU/Hr-ft. 2 ), 
maintain a constant temperature of -423°F at the inner 
tank wall, prevent frost build-up and air liquefaction on 
the surface of the Cl (exterior surface temperature greater 
than 32°F), and maintain temperatures above the 
minimum use temperatures of the various materials 
used in the tank wall, Cl, adhesives, and TPS. The 


only criteria for the re-entry/soak-through case were to 
not exceed the maximum, or fall below the minimum, 
use temperatures of the materials through-the-thickness. 
Temperatures may decrease for the re-entry/soak-through 
case due to aerothermal cooling. 

The results from both thermal analyses cases were 
evaluated to determine if the temperature limits of the 
materials were violated. If a temperature in a layer was 
near the maximum or minimum temperature limit of 
the material for either the ground-hold or re-entry/soak- 
through analyses, the insulation thicknesses of the Cl or 
TPS were decreased or increased accordingly and both 
analyses repeated. The thickness of the tank wall, 
adhesive layers, coatings, and blanket over wrap, were 
held constant. 

Insulation Concepts 

Two design concepts for insulation systems were 
analyzed in the thermal sizing study and are listed in 
Table 1. Temperature limits for these materials are 
shown in Figure 3, and the densities and temperature 
limits are listed in Table 2. Concept I had an A1 2219 
metallic alloy tank wall with a hybrid Cl bonded to the 
external surface with Lord 212™ adhesive. The hybrid 
Cl was a layer of PDL-1034 Cl and a layer of TEEK Cl 
[3]. In an actual RLV Cl, the foam of the hybrid Cl 
would be inside a honeycomb core to increase the 
tension strength, bondline strength, and robustness of 
the CL Unfortunately, thermal properties for the hybrid 
Cl in honeycomb core were not available for this study. 
The neat Cl thermal properties were used instead. 
Conformable reusable insulation (CRI) TPS was 
bonded to the hybrid Cl with RTV-560 adhesive. The 
CRI consisted of CRI batting insulation over-wrapped 
with the same fabric from an advanced flexible reusable 
insulation (AFRSI) TPS. The external surface or outer 
mold-line (OML) was thicker than the inner surface or 
inner mold-line (IML). The OML of the AFRSI fabric 
was coated with C-9 ceramic coating. Concept II had 
the same tank wall and Cl configuration as Concept I, 
but had a coated flexible reusable surface insulation 
(FRSI) TPS bonded to the external surface of the Cl 
with RTV-560. The OML of the FRSI was coated with 
hardened FRSI (HFRSI). Concept I was used on the 
windward or bottom surface, as shown in Figure 1, and 
Concept II was used on the leeward or upper surface. 

Thermal Load Cases 

Ground-hold and re-entry/soak-through load cases 
provided the temperature conditions to size the 
thickness of the Cl and the internal insulation of the 
TPS. The ground-hold condition case was developed to 
simulate the point in time when the vehicle is on the 
launch pad being filled with fuel just prior to launch. 
The ascent phase was not analyzed. The re-entry and 
soak-through conditions are combined because these 
conditions take place consecutively during the return to 
Earth phase of the vehicle. During re-entry, the exterior 
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surface of the vehicle experiences the highest 
aerothermal heating which was greater than the ascent 
phase. At soak-through, the TPS releases the heat 
energy that was absorbed during re-entry to the 
atmosphere and internal structure. Thus, the internal 
structure will possibly experience the maximum 
temperature after the vehicle has landed [4]. 

Table 1. Concepts, initial thicknesses, and upper 
and lower limits used in the thermal sizing study. 

Thickness 


(in) 


Concept 


Material 

Initial 

Lower 

limits, 

upper 

limits 



AFRSI OML 

0.027 

•> 


CRI* 

CRI batting 

1.5 

0.25, 2.0 



AFRSI IML 

0.009 

•> 

I 

RTV-560 f 

0.0075 



TEEK* 

0.75 

0.25, 1.5 


PDL-1034* 

0.75 

0.25, 1.5 


Lord 212™* 

0.0075 

•> 


Al 2219* 

0.16 

•> 


FRSI 

HFRSI f 

0.0075 

•> 


FRSI f 

1.50 

0.125, 1.5 

n 

RTV-560 f 

0.0075 


PDL-1034* 

0.75 

0.25, 1.5 


Lord 212™* 

0.0075 

•> 


Al 2219* 

0.16 

? 


Note: 

* Boeing proprietary data 
' Data is from [5] 


(T max s 1,800°F (CRI), 900°F (FRSI), 


mm. 


aO°F) 


Umax. * 450°F, 

Tmin. a -160°F) 

(Tmax. * 250°F, 
T min . a -460°F) 


(T max . s 250°F, 
Tmin. 2= -460°F) 



C-9 Ceramic, HFRSI 
coating 

CRI or FRSI TPS 
RTV-560 

TEEK Cl 1 50/50 
PDL-1034 Cl Hybrid Cl 


Lord 212™ 
Isogrid Al 2219 


2 

BTU/Hr.-ft. -°F [4] was used. Radiation to atmosphere 
was the boundary condition used for the re-entry/soak- 
through case with a heat flux load from the aerothermal 
loading. The normalized windward and leeward heat 
flux loads and pressure profiles for the re-entry/soak- 
through case are shown in Figure 5. The heat fluxes are 
normalized with respect to the maximum windward heat 
flux. 
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Figure 4. Two cases for the thermal sizing study: 
(a) ground-hold and (b) re-entry/soak-through. 


Table 2. Density and temperature limits of the 
material used in the thermal sizing study. 


Material 


Density Temperature (°F) 
(lbs./in. 3 ) Min. Max. 



AFRSI OML 

0.03484 

-200 

2,600 

CRI 

CRI batting 

0.006134 

-200 

1,800 


AFRSI IML 

0.03484 

-200 

2,600 

FRSI 

HFRSI 

0.05613 

-200 

1,160 

FRSI 

0.003125 

-200 

900 

RTV-560 

0.05093 

-160 

550 


TEEK 

0.001157 

-450 

450 

PDL-1034 

0.001447 

-450 

250 

Lord 212™ 

0.03935 

-450 

250 


Al 2219 

0.1030 

-450 

250 


ABAOUS Analysis 


Figure 3. Graphical depiction of the tank cross- 
section with minimum and maximum use 
temperature limits for the respective materials used 
in the thermal sizing study. 

The thermal boundary conditions in the model are 
depicted in Figure 4. A constant temperature of 
-423°F was maintained on the inner tank wall surface 
for the ground-hold case. An adiabatic thermal 
boundary condition was assumed at the inner tank wall 
surface due to the large volume of the tank for the re¬ 
entry/soak-through case. Natural convection to air was 
the boundary condition used on the exterior surface for 
the ground-hold case. A convection coefficient of 2.4 


A one-dimensional FEM model was used to model 
the through-the-thickness layers of the tank wall. A 
nonlinear transient heat transfer analysis was employed 
for the following studies: mesh convergence, ground- 
hold, and re-entry/soak-through. The model was 
analyzed in increments of 10 sec. The mesh 
convergence study was reported in a previous paper by 
Smeltzer et al. [6]. 

The initial wall designs had a total thickness of 3.2 
inches. Mesh patterns incorporating 261 nodes and 260 
elements for Concept I and 221 nodes and 220 elements 
for Concept II were used. The one-dimensional solid 
link element, DC1D2 in ABAQUS, was used in all the 
models analyzed. The elements were oriented in the 
through-the-thickness x-direction with a cross-sectional 
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area of 1.0 in. in the yz-plane. The boundary 
conditions varied as required by the particular analyses 
being performed. 



Time, Sec. 

Figure 5. The normalized heat fluxes on the 
windward and leeward centerlines of the vehicle 
and exterior atmospheric pressure during re-entry 
used in the thermal sizing study. 

The thermal conductivities for CRI and FRSI were 
both temperature and pressure dependent, however, the 
specific heats were only temperature dependent. 
Temperature and pressure dependence effected the 
thermal conductivity by reducing the thermal 
conductivity as the pressure and temperature decreased. 
The thermal conductivities and specific heats of the 
other materials, except for the coatings, were only 
temperature dependent. The coatings affected the value 
of the emissivities. The emissivities were held at 
constant values of 0.8 and 0.9 for Concepts I and II 
respectively due to constraints in ABAQUS, even 
though the emissivities were temperature dependent. 

The ground-hold case had an initial isothermal 
temperature of 50°F with an ambient air temperature of 
50°F to simulate a launch on a cold day. A 
temperature condition of -423°F was imposed on the 
inner tank wall. An initial isothermal temperature of 
100°F through-the-thickness was imposed on the 
models at the start of the re-entry/soak-through analysis 
case to simulate a warm initial temperature. An 
ambient temperature condition of 80°F was imposed on 
the exterior node of the models to simulate a 80°F 
temperature at the landing site. 

Results 

The thermal analysis and sizing results for both 
concepts of the ground-hold and re-entry/soak-through 
cases are presented in Figures 6 through 13, and Tables 
3 and 4. The variation of temperature results through- 
the-thickness of the tank wall for both concepts is 
plotted in the figures for selected time steps until 
steady-state temperatures for the ground-hold case were 
reached or a clear decrease in temperature for the re¬ 
entry/soak-through case began. The iterative results for 
the insulation thickness and areal weights are listed in 
Table 3. The through-the-thickness layers depicted in 
Figure 3 are indicated at the top of the plots in Figures 


6, 8, 10, and 12. Plots of the temperature results at the 
nodal interfaces at the top of each layer are shown in 
Figures 8, 10, 12, and 14. 

In Table 3, the final thickness results from 
ABAQUS are listed for Concepts I and II. The 
through-the-thickness temperatures are listed in Table 4. 
Eight iterations were required to obtain a feasible 
solution for Concept I and six iterations were required 
for Concept II. The final areal weights for the 
insulation are 0.01112 lbs./in. 2 and 0.004093 lbs./in. 2 
for Concepts I and II respectively. 

Table 3. Cl and TPS insulation thickness and 


total weight results. 

Concept 

Iter. 

no. 

CI 

(in.) 

TPS 

(in.) 

Areal wt. T 
(lbs./in. 2 ) 


1 

1.5 

1.5 

0.01309 


2 

1.2 

1.2 

0.01086 


3 

1.0 

1.4 

0.01182 

I 

4 

0.8 

1.4 

0.01156 


5 

0.8 

1.3 

0.01095 


6 

0.7 

1.3 

0.01143 


7 

0.7 

1.35 

0.01112 


1 

1.5 

1.5 

0.007739 

n 

2 

1.8 

0.25 

0.004223 

3 

1.6 

0.125 

0.003963 


4 

1.7 

0.125 

0.004093 


Note: 


Cl thickness is evenly divided between the PDL-1034 and TEEK. 

TPS thickness is the CRI batting or FRSI for Concepts I and II 
respectively. 

f The areal weights for Concept I includes the areal weights of the 
Lord 212™, Cl, RTV-560, AFRSIIML, CRI batting, AFRSI OML, 
and C-9 coating. The areal weights for Concept II includes the areal 
weights of the Lord 212™, Cl, RTV-560, FRSI, and HFRSI. 

The ground-hold results plotted in Figures 6 
through 8 indicate that both concepts will operate 
efficiently and effectively with 0.7-in.-thick and 1.7-in.- 
thick Cl for Concepts I and II, respectively. The results 
in Figures 6 and 8 display how the concepts cooled 
over time with a -423°F tank wall inner surface. A 
steady-state temperature condition was achieved for 
Concept I by t = 6,000 sec. in Figure 7, with a surface 
temperature of 38.18°F. Concept II reached a steady- 
state temperature condition by t = 5,000 sec. in Figure 
9, and had a surface temperature of 35.42°F. Both 
concepts maintained a temperature above -160°F at the 
RTV-560 bondline and the upper surface should not 
develop frost (the exterior temperature is greater than 
32°F** in Figures 6 through 9). The maximum average 
through-the-thickness heat flux was 28.4 BTU/Hr.-ft. 2 
and 33.9 BTU/Hr.-ft. 2 for Concepts I and II, 


* * A surface temperature of 32°F does not ensure that frost build¬ 
up will not occur. Other factors such as humidity, ambient air 
temperature, radiation, and atmospheric pressure dictate whether 

frost will form on a surface. The boundary conditions and initial 
temperatures imposed in this study are conservative and will 
minimize the potential of frost formation on the actual surface. 
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respectively. Both concepts had an average heat flux 
less than the maximum limit of 45 BTU/Hr.-ft. 2 to 
prevent LH 2 boil-off. 



Depth, in. 

Figure 6. Through-the-thickness temperature 
results for Concept I with CRI TPS for the 
ground-hold case. 



Figure 7. Temperature results at the upper surface 
of each layer of Concept I with CRI TPS for the 
ground-hold case. 


The re-entry/soak-through temperature results from 
ABAQUS are plotted in Figures 10 through and 13. 
The propagation of the temperature pulse through-the- 
thickness is shown in Figures 10 and 12 for Concepts I 
and II, respectively. The results in Figure 11 indicate 
that Concept I achieved maximum temperature 
conditions at the critical layers (tank wall and Cl to 
TPS interface) at approximately 2,100 sec. to 2,300 
sec. The results in Figure 13 indicate that Concept II 
achieved the maximum temperature conditions at the 
same critical layers at approximately 1,200 sec. to 
1,400 sec. The maximum temperature on the exterior 
surface was 1,676°F at 1,200 sec. for Concepts I and 
542.6°F at 1,200 sec. for Concepts II. The through- 
the-thickness temperatures began to decrease in Figure 
11 at 2,300 sec. for Concept I and in Figure 13 at 
2,900 sec. for Concept II. The through-the-thickness 
temperatures, listed in Table 4, were below the 
maximum use temperatures of the materials shown in 
Figure 3 for both concepts. There was also a 


significant decrease in temperature on the upper nodes 
for both concepts in Figures 10 and 12. This was due 
to simulated aerothermal cooling as the vehicle passed 
through the atmosphere. 


AI2219 PDL-1034 TEEK FRSI 



Figure 8. Through-the-thickness temperature 
results for Concept II with FRSI TPS for the 
ground-hold case. 



Figure 9. Temperature results at the upper surface 
of each layer of Concept II with FRSI TPS for the 
ground-hold case. 


Concept I required a greater number of iterations 
than Concept II to achieve an optimal design. The 
additional iterations were precipitated from the Cl and 
TPS working as a system during the ground-hold and 
re-entry/soak-through analyses. The TPS acted as a 
cryogenic insulator during ground-hold and the Cl acted 
as a high temperature insulator during re-entry/soak- 
through. As the insulation thickness was reduced, low 
and high temperature limits were violated, and the 
insulation thicknesses had to be revised. The minimum 
temperature at the RTV-560 bondline and prevention of 
frost on the exterior surface dominated the sizing of the 
insulation for the ground-hold case for Concept I. The 
prevention of frost and the maximum average heat flux 
through-the-thickness for the ground-hold case dictated 
the thickness of the Cl for Concept II while the FRSI 
was driven to the minimum thickness allowed. Since 
the FRSI was at its minimum thickness, the FRSI did 
not have a large contribution as a cryogenic insulator, 
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which required Concept II to have almost twice as 
much Cl as Concept I. 



Figure 10. Through-the-thickness temperature 
results for Concept I with CRI TPS for the re¬ 
entry/soak-through case. 



Figure 11. Temperature results at the upper surface 
of each layer of Concept I with CRI TPS for the re¬ 
entry/soak-through case. 


STRUCTURAL SIZING STUDY 

The skin thicknesses and isogrid geometry of the 
tank wall were sized for an aluminum tank using the 
pressure and mechanical loads from two conditions: 
proof-pressure testing and ground handling. Before a 
tank is brought into service, proof-pressure testing will 
be performed at 1.05 times the maximum pressure. 
Only the barrel section of the cryogenic tank was 
analyzed in this study to simplify the analyses. 

Sizing Approach 

The design and optimization of metallic stiffened 
shells for buckling and stress constraints is 
characterized by a large number of local optima. 
Because the design space is characterized by multiple 
local optima, global optimization plays an important 
role. A deterministic optimization scheme was 
developed using the divide rectangles (DIRECT) [7, 8] 
global optimization algorithm in conjunction with 
PANDA2 [9], generating a number of feasible or nearly 


feasible designs for reliability analysis using Monte 
Carlo simulation (MCS). 

AI2219 PDL-1034 TEEK FRSI 



Figure 12. Through-the-thickness temperature 
results for Concept II with FRSI TPS for the re¬ 
entry/soak-through case. 


Table 4. Results for minimum and maximum 
temperatures at the upper surface of each layer. 


Concept 

Layer 

Ground- 

hold 

(F°) 

Re-entry/ 

Soak- 

through 

(F°) 


AFRSI OMF 

38.2 

1,676 


CRI batting 

34.7 

1,663 


AFRSI IMF 

-153.8 

345.3 

I 

RTV-560 

-154.1 

344.5 

TEEK 

-154.2 

344.1 


PDF-1034 

-240.6 

246.4 


Lord 212™ 

-422.9 

177.0 


A1 2219 

-423.0 

177.0 


HFRSI 

35.42 

542.6 


FRSI 

35.15 

539.6 


RTV-560 

18.41 

428.4 

n 

TEEK 

18.32 

427.7 


PDF-1034 

-162.5 

249.1 


Lord 212™ 

-422.8 

141.1 


A1 2219 

-423.0 

141.1 


PANDA2 Program 

PANDA2 is an analysis and optimization program 
developed by David Bushnell [9] that is capable of 
analyzing and optimizing flat or radially symmetric 
panels stiffened with combinations of internal and 
external stringers, ring frames, and grids of various 
cross-sections. PANDA2 enables the analysis by using 
a combination of relatively simple models to capture 
individual failure modes. Possible failure modes 
include: yielding, and general stiffener and skin pocket 
buckling. Optimization is performed using gradient- 
based algorithms [10], with mechanical failure 
constraints and user-defined geometric constraints. 
Feasibility is determined using constraint margins, with 


6 

American Institute of Aeronautics and Astronautics 























































































AIAA Paper 2004-1931 


the convention that negative margins are infeasible. 
Mechanical failure constraints are defined such that 
maximum violation corresponds to a margin of -1. 
Geometric margins are evaluated using a user-defined 
inequality, and may in fact be much lower than -1 in 
various regions of design space. 



explanation of this procedure is given in references [7] 
and [8]. 


w 2 = 358.2 



O 

67207 _ 


o 

O 

O 

352.4 

. 600_4_ 

O 

358.2 

200.6 


200.6 


- 1 st division 

— — 2nd division 



Potentially 
optimal box 









Figure 14. (a) 1 st iteration of DIRECT which shows 
windowing for a 2-dimensional example, the 
Goldstein-Price function and (b) further iterations. 


Figure 13. Temperature results at the upper surface 
of each layer of Concept II with FRSI TPS for the 
re-entry/soak-through case. 

PANDA2 provides a single deterministic optimum 
design, and cannot directly calculate the probability of 
failure (POF) of a given design, making reliability 
optimization difficult. A new optimization algorithm, 
DIRECT [9, 10], was introduced in order to facilitate 
reliability optimization, broadening the range of design 
problems to which PANDA2 may be applied. In 
addition, MCS and response surface approximation [11] 
were used for efficient reliability analysis. 

The critical margins (Mi, M 2 ) are defined as the 
most negative margin (that was associated with the 
most violated constraint) evaluated at each position on 
the shell. The critical margin, Mi, is evaluated either 
between the ring frames, if the ring frames are present, 
or else at the middle of the panel. The critical margin 
at the ring frames is M 2 . If no ring frames are present, 
M 2 is undefined and was assumed to be equal to zero. 

DIRECT Algorithm 

DIRECT, a global sampling optimization 
algorithm, was selected over gradient-based 
optimization because the designs space has many local 
optima. Additionally, changing critical constraints and 
imperfection-induced discontinuities caused sharp 
variations in the critical margin surface. 

DIRECT is a modified Lipschitzian optimization 
algorithm [7, 8] that begins by evaluating the function 
to be minimized at the center of the design space. The 
algorithm then divides the space into thirds based on 
evaluations of the objective function, and chooses 
potentially optimum boxes for further division based 
on the objective function and the size of the box. This 
division operation is repeated, as illustrated in Figure 
14. Given infinite iterations, DIRECT will evaluate 
every point in the design space. A more detailed 


Implementation 

A penalty function was used for infeasible designs. 
A more detailed description of the implementation of 
the penalty function is in reference [12]. Ideally, this 
penalty, W pen , is the increase in weight, W , required to 
attain feasibility. To eliminate designs with large 
negative margins, a penalty function with a quadratic 
term was used. The penalty function also contains a 
linear term to force the optimum to favor the feasible 
domain over slightly infeasible points. Thus, a suitable 
penalized weight function is of the form 

w pe „ = w[i + + AM,) + a(« 2 m 2 + AM 2 )] (!) 

where a and |3 are penalty parameters, and 6i is a 
switch, equal to 0 when Mi is greater than 0 (feasible) 
and 1 when Mi is less than 0 (infeasible). 

To simplify the problem, the penalty parameters a 
and (3 were assumed to be, ai = a 2 = a and |3i = |3 2 = 
|3. It is sufficient and rather simple to find a and |3 by 
trial and error such that the best designs are either 
feasible or nearly so. Reasonable values for the types of 
problems used in this study are a = 4 and (3 = 1. 

Reliability Analysis 

MCS was used in the reliability analysis because 
this simulation technique is easy to implement, robust, 
and accurate with sufficiently large samples. However, 
MCS requires a large number of analyses to obtain an 
accurate estimate of small failure probabilities. 
Response surface approximations were used to reduce 
the number analyses, by fitting a closed-form 
approximation to the limit state function. Response 
surface approximations usually fit low order 
polynomials to the structural response in terms of 
random variables. Qu et al. [13] developed a global 
analysis response surface (ARS) approach in the unified 
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space of design and random variables to substantially 
reduce the number of response surface approximations. 

Aluminum Tank Optimization 

The tank barrel optimized was an aluminum 
cylinder (A1 2219), stiffened externally with J-shaped 
ring frame stiffeners, and internally with a blade-shaped 
isogrid oriented circumferentially as shown in Figure 
15. The material properties for A1 2219 are listed in 
Table 5. The length, Li, of the tank barrel was 300 in., 
and the radius, r, was 160 in. Using symmetry, it 
sufficed to model half of a cylindrical tank, spanning 
180°. The tank was subjected to both internal pressure 
and axial compression. The pressure partially relieved 
the influence of buckling due to axial compression, an 
effect known as pressure stabilization. The first load 
set was a proof-pressure of 35 psi with no axial load. 
The second load set was due to ground handling with 
an operating pressure of 5 psi and an axial compression 
load of 1,000 lbs./in. Thermal loads were not 
considered since the tank was assumed to be at 
isothermal conditions and the tank material was 
assumed to be isotropic. 


Table 5. A1 2219 material properties and 
coefficients of variation. 


Var. 

Density 

(p) 

(lbs./ft. 3 ) 

Young’s 

Modulus 

(E) 

(psi) 

Poisson’s 

ratio 

04 

Allow. 

Stress 

(Oa) 

(psi) 

Mean 

0.1030 

0.107x10 s 

0.34 

0.58xl0 5 

COY 

— 

0.03 

0.03 

0.05 


The optimization scheme varied the isogrid spacing, 
the isogrid blade height, and the thicknesses of the skin 
and isogrid blades. The ring frame geometry was 
considered fixed and listed in Table 6. The 
optimization was subject to PANDA2 failure 
constraints with the objective of minimizing the 
structural weight of the tank. 

The PANDA2 analysis was based on a procedure 
developed by Lamberti et al. [14], for the optimization 
of a similar aluminum tank with the ring frame 
geometry removed. The lower and upper bounds for 
optimization are given in Table 7 along with a design 
obtained by the internal optimizer in PANDA2. 


Table 6. Fixed tank geometry for the ring frame. 


Variable 

Fixed geometry 
(in.) 

Spacing (b ring ) 

113.9 

Base (bi) 

0.4 

Web height Ow) 

2.0 

Flange (w) 

0.48 

Base thickness (T 3 ) 

0.03 

Web thickness (T 4 ) 

0.06212 

Flange thickness (T 5 ) 

0.07952 


Table 7. Variable tank geometry with bounds and 
PANDA2 optimum design. 


Variable 

Lower 

bound 

(in.) 

Upper 

bound 

(in.) 

PANDA2 

results 

(in.) 


Spacing (b,, rld ) 

5 

15 

10.71 

Isogrid 

Height (h erid ) 

1 

5 

2.049 


Thickness (T 2 ) 

0.08 

0.25 

0.0989 

Skin thickness (Ti) 

0.03 

0.25 

0.1261 


Weight 

— 

— 

5,196 lbs. 


DIRECT Optimization 

A feasible optimum design was obtained using 
DIRECT with a = 4 and (3=1 and is listed in Table 8 
with results for other penalty parameters. The range of 
geometries for the designs with a penalized weight 
similar to the optimum (near-optima) are given in Table 
9 for the case of a = 4 and (3 = 1. These designs are 
fairly wide-spread around the design space. 
Approximately 1,000 function evaluations were needed 
to converge to the design in Table 10 that was 4 percent 
lighter than the PANDA2 optimum. DIRECT 
optimization did not evaluate any designs similar to the 
optimum determined in PANDA2, although additional 
iterations of DIRECT would eventually reveal designs 
similar to the PANDA2 optimum. 



Figure 15. The basic tank barrel geometry for an 
isogrid stiffened aluminum cryogenic tank with 
exterior ring frames. 

Accurate critical margin response surfaces were 
constructed to evaluate the reliability of a design point. 
The uncertainty in the material properties was 
represented by three random variables, two for elastic 
properties (modulus of elasticity, E, and Poisson’s 
ratio, fi) and one for the allowable stress (c^), which 
were assumed to be normally distributed and 
uncorrelated. The mean values and coefficients of 
variation of the uncertainties in material properties are 
shown in Table 5. The four design variables ( b , h , 77 
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and T 2 ) also have randomness in them due to 
manufacturing uncertainties and were assumed to be 
uniformly distributed. Table 11 shows the percentage 
variation of the random design variables. The data in 
Table 11 are obtained based on rough estimates, and are 
intended for illustration purposes only. 


Table 8. DIRECT optimum designs with various 
penalty coefficients. _ 


Penalty 

Penalized 
Weight weight 

(lbs.) (full tank) 

Mi M 2 

a |3 

(lbs.) 


1 0 

4,764 4,821.6 

0.197 -0.110 

2 0 

4,724 4,813.1 

0.179 -0.097 

3 0 

4,760 4,834.6 

0.190 -0.072 

4 0 

4,944 4,945.1 

0.245 -0.007 

4 1 

4,970 4,970.0 

0.177 0 

Table 9. Range of design variables for designs with 
Wpen within 200 lbs. of optimum full tank weight. 

Entity ^ ^ 

Tr T 2 

(in.) (in.) 

Low 

5.021 1.617 

0.0936 0.1197 

High 

6.187 2.012 

0.1104 0.1291 

% of bounds 11.7% 9.9% 

9.9% 4.3% 

Table 10. Comparison of DIRECT (a = 4, |3 = 1) 
and PANDA2 optima. 


Method 

Variable DIRECT 

optimum 

PANDA2 

optimum 

b 

5.043 in. 

10.71 in. 

h 

1.617 in. 

2.049 in. 

T, 

0.0957 in. 

0.0986 in. 

t 2 

0.1197 in. 

0.1261 in. 

Weight 

4,970 lbs. 

5,196 lbs. 

Mi 

0.239 

-9.95xl0’ 6 

m 2 

0.0 

0.0189 

Table 11. Percentage variation of the random 
design variables around design value. 


b h 

Tr T 2 

% Variation ±1% ±2% 

±4% ±4% 


Response surface approximation approaches were 
employed to fit response surfaces to the most critical 


safety margins for the designs centered around the near¬ 
optima selected based on Latin Hypercube sampling and 
the variabilities. As in the deterministic optimization, 
isogrid panels were analyzed using PANDA2 with 
constraint safety factors adjusted to 1.0. The most 
critical margins of the two load cases were then 
extracted. Polynomial response surface approximations 
were fitted to the worst safety margins in terms of 
random variables. POF for the panel system and each 
load case were calculated by performing MCS on the 
polynomial response surfaces [13]. This calculation 
was performed for the PANDA2 and DIRECT optima 
given in Table 10, with response surface statistics given 
in Table 12, and POF given in Table 13. Both points 
were modeled with acceptable accuracy. Further 
investigation of the design space for feasible solutions 
was not necessary since the DIRECT design meets the 
very strict POF constraint, at most 1 in 10 8 . 


Table 13. POF calculated by MCS. 


Method 

No. of 
samples 
in MCS 

POF at 
mid-ring 
frame 

POF at the 
ring frame 

System 

POF 

DIRECT 

optimum 

lxlO 9 

0 

0 

0 

PANDA2 

optimum 

10xl0 6 

0 

394x1 O' 6 

394x1 O' 6 


COMBINED RESULTS 

The results from the thermal sizing study and 
structural study are combined in Table 14. An efficient 
and minimum areal weight for a barrel section is shown 
in Table 14 for the windward (Concept I) and leeward 
(Concept II) sides of the vehicle. The windward side of 
the vehicle is much heavier (34.2 percent) than the 
leeward side due to the additional TPS required to 
withstand the high heat flux load on the windward side 
during re-entry. The additional TPS, however, reduced 
the amount of Cl, since the TPS acts as a cryogenic 
insulator. 

To further minimize the weight for the whole barrel 
section, similar thermal and structural sizing studies 
would be required around the circumference of the 
barrel. Plug thermal sizing around the circumference 
would enable the tailoring of TPS from the windward 
side to the leeward side of the vehicle. Additionally, 
the tank structure would be separated into four sections 
such as the crown, two side panels, and keel in a more 
detailed structural sizing study. 


Table 12. ARS approximation to the worst margin using Latin Hypercube sampling. 


A • A 

Cubic ARS 


Quadratic ARS 

Quantity 

Mi fit 

M 2 fit 

Mi fit 

M 2 fit 

R 2 adi. 

0.9999 

0.9377 

0.9085 

0.9998 

RMSE predictor 

0.000253 

0.0416 

0.0261 

0.000603 

Mean of response 

0.6929 

1.2812 

0.4601 

0.1251 
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Table 14. Combined areal weight results from the thermal and structural sizing studies. 


Concept 

Tank 

weight from 
DIRECT 
(lbs./in. 2 ) 

Cl 

weight* 

(lbs./in. 2 ) 

TPS 
weight 
(lbs./in. 2 ) 

Total 
insulation 
weight 
(lbs./in. 2 ) 

Combined 
weights 
(lbs./in. 2 ) 

I 

0.01648 

0.001207 

0.009917 

0.01112 

0.02760 

n 

0.01648 

0.002509 

0.001584 

0.004093 

0.02057 


Note: 

* Cl areal weight includes the areal weights of the Lord 212™ adhesive and CL 

f TPS areal weights for Concept I includes the areal weights of the RTV-560, AFRSI IML, CRI batting, AFRSI OML, and C-9 coating. TPS 
areal weights for Concept II includes the areal weights of the RTV-560, FRSI, and HFRSI. 


CONCLUDING REMARKS 

Thermal and structural analysis tools were 
developed for the sizing of insulation and structural 
geometry of integrated cryogenic propellant tank 
concepts for a reusable launch vehicle (RLV). These 
tools were used to size the thickness of the cryogenic 
insulation (Cl), insulation thickness for thermal 
protection system (TPS), and isogrid stiffener geometry 
and skin thickness of an aluminum 2219-T87 tank 
wall. The analytical studies and results described in 
this paper are a portion of the basic information 
required for the development of light-weight reusable 
integrated cryogenic propellant tanks. 

A thermal sizing study to determine the thickness 
of the insulation was successfully performed using 
ABAQUS, a commercial nonlinear finite element 
method (FEM) analysis code for separate but related 
analytical studies for ground-hold and re-entry/soak- 
through conditions. The nonlinear transient heat 
transfer analyses used the temperature and pressure 
dependent thermal material properties to determine the 
through-the-thickness temperatures of two concepts. 
Concept I employed conformable reusable insulation 
(CRI) TPS for the windward side while Concept II, the 
leeward side, used flexible reusable surface insulation 
(FRSI) TPS. Concept I used approximately half of the 
Cl that Concept II required but weighed 34.2 percent 
more in combined insulation weight. The increased 
weight was due to the additional TPS required to resist 
re-entry aerothermal heating, while Concept II had the 
thinnest layer of FRSI possible. The prevention of 
frost accumulation on the external surface for the 
ground-hold case dominated the sizing for both 
concepts. 

Alternative structural designs were obtained by 
using the DIRECT optimizer to explore the design 
space, allowing for more accurate reliability 
optimization. The design obtained from DIRECT was 
4.3 percent lighter than the PANDA2 deterministic 
design and had a near zero probability of failure 
compared to the 3 64x1 O' 6 probability of failure for the 
PANDA2 design. 

The methods and analytical tools developed for 
these studies provided an efficient means to analyze, 
size, and evaluate candidate-integrated designs for future 


RLVs. The thermal sizing study established a method 
for sizing the thickness of insulation layers using FEM. 
The structural sizing study developed an optimization 
method that produced more robust alternate structural 
designs than deterministic optimization methods by 
using reliability optimization techniques in conjunction 
with deterministic optimization techniques. 
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